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Abstract 

In this paper we address the questions of perfectly sampling a Gibbs measure with 
infinite range interactions and of perfectly sampling the measure together with its finite 
range approximations. We solve these questions by introducing a perfect simulation 
algorithm for the measure and for the coupled measures. The algorithm works for 
general Gibbsian interaction under requirements on the tails of the interaction. As a 
consequence we obtain an upper bound for the error we make when sampling from a 
finite range approximation instead of the true infinite range measure. 
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1 Introduction 

Consider a Gibbs measure with an infinite range potential and its finite range approxima- 
tion obtained by truncating the interaction at a certain range. Can we construct a perfect 
sampling of the Gibbs measure with an infinite range potential? Further, if we make a 
local inspection of a perfect sampling of the finite range approximation, how often does it 
coincide with a sample from the original infinite range measure? We address these ques- 
tions by introducing a new perfect simulation algorithm for the original Gibbs measure 
and a coupled perfect simulation algorithm for the original measure together with its finite 
range approximation. 

By a perfect simulation algorithm we mean a function mapping a family of independent 
uniform random variables on the set of configurations having as law the infinite range 
Gibbs measure. By a coupled perfect simulation algorithm we mean a function mapping 
a family of independent uniform random variables on the set of pairs of configurations 
having as first marginal the original infinite range measure and as second marginal the 
finite range approximation. These functions are defined through two coupled Glauber 
dynamics having as reversible measures the original Gibbs measure and its finite range 
approximation, respectively. The perfect simulation algorithm produces coupled samples 
of the two measures. 

The main difficulty in this constructive approach is that we have to deal with infinite range 
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spin flip rates. This difficulty is overcome by a decomposition of the spin Hip rates of the 
Glauber dynamics as a convex sum of local spin flip rates. 

Our results are obtained under high temperature conditions and assuming that the inter- 
action decays fast enough. We do not need any monotonicity property of the infinite range 
Gibbs measure. Under these conditions, we show that the algorithm that we propose stops 
almost surely, after a finite number of steps. Moreover, we give an upper bound, uniform 
in space, for the probability of local discrepancy between the coupled realizations of the 
two measures. This is the content of Theorem [1] As a corollary. Theorem U] shows that 
the same upper bound holds for Ornstein's d-distance between the two measures. 

An upper bound for the d-distance slightly worse than the one obtained in Theorem H] 
can be obtained under less restrictive assumptions, by applying Dobrushin's contraction 
method. We refer the reader to the series of papers of Dobrushin (1968-1970) as well as 
to Presutti (2009) for a recent nice self-contained presentation of the subject. However, 
the contraction method is not constructive and does not provide an explicit sampling 
procedure for the measures. 

The approach followed in the present paper was suggested by a recent work of Galves et 
al. (2008). There, the ffip rates of an infinite range interacting multicolor system were 
decomposed as convex combinations of local range ffip rates. With respect to this work, 
the novelty of what follows is that the decomposition is explicit and more suitable in view 
of the intended coupling result. The idea of decomposing ffip rates can be traced back 
to Ferrari (1990) in which something with the same flavor was done for a finite range 
interacting particle system and to Ferrari et al. (2000) in which the probability transition 
of a stochastic chain of infinite order was decomposed as a convex combination of finite 
order probability transitions. We refer to Galves et al. (2008) for more details and related 
literature. 

This paper is organized as follows. In Section 2, we present the basic notations and the 
main result. Theorem [H In Section 3, we introduce the Glauber dynamics and the asso- 
ciated processes. The representation of the spin flip rates as convex combination of local 
rates is given in Theorem [2j In Section 4.1 we introduce the backward sketch process 
which is the basis of the perfect simulation algorithm and collect the main properties of 
this process in Proposition [3l Section 4.2 presents the coupled perfect simulation algorithm 
for the couple of measures. In Section 5 we prove Theorem [TJ We conclude the paper 
in Section 6 discussing how certain features of high temperature Gibbs systems can be 
obtained applying our construction. In particular we give results on Ornstein's d-distance 
between the infinite range and the finite range Gibbs measure, and on the decay of corre- 
lations. These results are obtained under requirements stronger than the one needed using 
the Dobrushin contraction method, since they are based on our explicit coupled perfect 
simulation procedure. 
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2 Definitions and main results 

Let A = {—1, 1} and S = A^"^ be the set of spin configurations. We endow S with the 
product sigma algebra. We define on Z'^ the norm, ||i|| = Ylk=i l^^l- 

Configurations will be denoted by Greek letters a,r], A point i e will be called a 

site. For any i G Z**, notations a{i) and ai will be used in order to denote the value of 
the configuration a at site i. By extension, for any subset V C Z,'^, ^"(y) and ay € A^ 
denote the restriction of the configuration a to the set of positions in V. If V C Z^ is a 
finite subset, we denote V (£ U^. In this case, we write \y\ for the cardinal of V. Finally, 
for any z G Z*^ and a E S, we shall denote cr' the modified configuration 

= cr{j), for aU j ^ i, and a\i) = -a{i). 

Definition 1 An interaction is a collection {Jb,B Z*^} of real numbers indexed by 
finite subsets B ^ U^, \B\ > 2, which satisfies 

sup J2 \M < oo. (2.1) 

i^^'' B-.ieB 

For finite B C Z'^ and a e S let 

XB{cr) = JJ(j(i). 

Definition 2 A probability measure on S is said to be a Gibbs measure relative to the 
interaction {Jb,B (E Z'^} if for all i G Z*^ and for any fixed ( E S, a version of the 
conditional probability /x({(7 : a{i) = C(,i)Wij) = C(i) for cdl j 7^ i}) is given by 

^i{{a : a{i) = C(i)|a(j) = C(j) for all j / z}) = — — (2.2) 

1 + exp(-2 }ZB:ieB JbXb{Q) 

To define the Gibbs measure /xt^' relative to the interaction truncated at range L we set 

Bi{L) = {j eZ'^:\\j-i\\<L}. (2.3) 

Then /i'^l is the probability measure on 5 such that for all i £ Z^ and for any fixed C ^ 
a version of the conditional probability /x[^1({(T : a{i) = Ci'i-)\(^{j) = C(j) foi^ j ¥^ i}) is 
given by 

fil'^^ia : a{i) = C{i)\a{j) = ({j) for all j + i}) 

1 

1 + exp(-2^^^.gB,BcBi(L) JbXb{0) ' 

We consider the interaction = P Jb where /? is a positive parameter. The associated 
Gibbs measures will be denoted by /x and omitting to write the explicit dependence 
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on 13. We denote by {X{i),i E Z'^), € Z'^) two random configurations on S 

distributed according to /i and respectively. 

All processes that we consider in this paper will be constructed as functions of an inde- 
pendent family of uniform random variables ([/„(i),n G N, i G TJ^). Denote {Vl,A^P) the 
probability space on which is defined the independent family of uniform random variables. 

Definition 3 A perfect sampling algorithm of the measure fi is a map F from [0, l]^^^'* 
to {-1,+!}^'*, such that 

F{Un{j),n E N,j E 'L'^) has distribution fi. 

A coupled perfect sampling algorithm for the pair {{X{i),X^^^{i)),i E Z*^) is a map G from 
[0,1]^^^" to such that 

G{Un{j),n E N,i E Z'^) has marginals {X{i),i E Z'') and {X^^\i),i E Z'^'). 

Since F produces an unbiased sample from the target distribution the algorithm is 
denoted "perfect". 

Definition 4 The perfect sampling algorithm F of the measure fi is said to stop almost 
surely after a finite number of steps if for every site i E Z^, there exists a finite subset F^^^ 
of Z and a finite random variable X^j.Qp > 1 such that if 

Kij) = UnU) for all j E , 1 < n < iV« ^p, (2.4) 

then 

F{Un{j),n E N,i E Z'^)(i) = F(C/;(j),n E N,i E Z'^)(i). (2.5) 

Note that (p3]) implies that F{Unij),n E N, j E Z'^)(i) depends only on C/„(j) for 1 < n < 
^STOP- '^^^ following theorem is our main result. 

Theorem 1 Assume that 

supJ^|A(A;)I |Jb||<oo. (2.6) 

1. There exists /3c > such that for any (3 < fie there is a perfect sampling algorithm 
of ji. This algorithm stops after a finite number of steps almost surely and 

sup E{N^%Op)<-^ (2-7) 

ieZ'i 7 

where 7 is given in (j4.37p . 

2. Moreover, for any (3 < f3c, there exists a coupled perfect sampling algorithm for the 
pair {{X{i),X^^\i)),i E Z"^) satisfying 

sup P(X(i) / X[^l(i)) < - sup (1 - e~'^^B:.eB,B^B,iL) \Jb\\ _ (3.8) 
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Remark 1 The constant Pc the solution of equation (|5.39|) in Section [3 



Remark 2 In case of an interaction defined through a pairwise potential {J{i,j), {i,j) G 
^ ^a!| ^^^^ J{i,i) = for i G Z"', condition (|2.6p of Theorem\^ reads as 

sup^|i?.(A;)| <°°- 

k \3:¥-,\\=k j 

In this case, 

supP(X(i) ^ X^^\i)) < - sup fl - e"'^^^^ll-^ll>^l-^(^'^')l) . 

The above perfect sampling algorithms are based on a coupled construction of two pro- 
cesses having /i and yut^^ as reversible measures. These two processes will be introduced 
in the next section. The definition of the algorithms can only be given in Section [H The 
proof of Theorem [T] is given in Section [SJ 



3 Glauber dynamics 

We now introduce a Glauber dynamics having /x as reversible measure. This is an inter- 
acting particle system {at{i),i G Z'^jt G M) taking values in S. Sometimes we will also use 
the short notation (o"j)t for the interacting particle system. To describe the process, we 
need some extra notation. 

For any i G Z*^, a G {—1, 1}^'*, and /3 > 0, we define Ci{a), which is the rate at which the 
spin i flips when the system is in the configuration a, as 

Ci{a) = exp I -/? ^ JBXBicr) | • (3.9) 
\ BA&B / 

The generator Q of the process {o't)t is defined on cylinder functions / : 5" ^ M as follows 

gf{a) = Y,Ci{a)[f{a')-f{a)]. (3.10) 

Observe that the condition (j2.6p of Theorem [1] implies that 

sup ^ \B\ \Jb\ < oo, 

* B:ieB 

where \B\ designs the cardinal of the set B. Hence the rates Ci{a) are uniformly bounded 
in i and a, and 

sup sup \ci{a) — Ci{a^)\ < oo. (3-11) 

Therefore, Theorem 3.9 of Chapter 1 of Liggett (1985) implies that G is the generator of a 
Markov process (o"t)t on SS. By construction, the process {crt)t is reversible with respect 
to the Gibbs measure fi corresponding to the interaction = (3 Jb- 
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For any ^ > 1, we now define a Glauber dynamics having /xl^l as reversible measure. More 

e"i fW) = E 4'V)[/(<^") - K")], (3.12) 



precisely, we consider the Markov process [a^ )t on SS having generator 



where the rates cf^ are given by 



cf{a)=rfe^Y>\-(i E JBXB{a)\ (3.13) 

and 

Let us stress the fact that compared to the usual definition of a Glauber dynamics, an 
extra factor rf^ appears in the definition of the rates cf^ . This extra factor is important in 

view of the intended coupling of both processes at and af^ (see Theorem [2] below) and does 
not change the equilibrium behavior of the process as is shown in the next proposition. 

Proposition 1 The process {<J^)t is reversible with respect to the Gibbs measure fi^^^ 
relative to the pairwise interaction Jg truncated at range i. 



Proof By construction, we have that 



cfl(a) _/.M(a^ 



c:- [a" 



Thus, by Proposition 2.7 of Chapter IV of Liggett (1985), it follows that (crf')t is reversible 
with respect to //t^^. • 

The main tool to prove Theorem [1] is a coupled construction of the processes {at{i),i G 
Z"^, t G M) and (erf (i), i G Z"^, t G M) for ^ = L, using the basic family of uniform random 
variables {Un{i),n G N, i G Z'^). This coupling is based on a decomposition of the spin flip 
rates Ci{a) as a convex combination of local range spin flip rates. 

To present this decomposition we need the following notation. We define 

Mi = 2e'^'^B:ieB\JB\^ 

and the sequence {\i{k))k>o as follows 

e-'^(^T,B■.^eB\JB\ if fc = 0, 

Aj(A;) = < e~^^s^'-e-s.sz:si(i) '-^-bI _ g-2/3Es:ies \Jb\ if k = 1, 

Q~P^B:i£B,B(tBi{k) l-^sl _ ^—P^B:i£B,B(tBi{k-l) I'^sl if A; > 1 

Finally, for all cr G S" and for any k > 2, we define the update probabilities 



J''^(_Cr(i)\a) = —e'^^B■.^SB,BCBi(k-l)JBXB{<T) 



1 

,BCBi(k),Bg^Bi(k-l) JBXB{f^) _ ^-I^T,B:iSB,BCBi(k),B(^Bi{k-l) \-^b\ 
\ _ Q~^^B:ieB,BCBi{k),B^Bi{k-l) l-^sl 
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and for k = 1, 



To obtain a probability measure on A, we define for all A; > 1, 

pf\a{i)\a) = l-pf\-a{i)\a). 

Finally, for A; = 0, we define 

= T^\-l) = l- (3.14) 

Since this last probability does not depend on the site i we omit the subscript i. Note 
that, by construction, for any integer A; > 1, the probabilities pf^\a\a) depend only on 
a{Bi{k)). 

Now we state the decomposition theorem. 

Theorem 2 Let us assume that the uniform summability condition Ii2.1\) holds. 

1. The sequence {Xi{k))k>o defines a probability distribution on the set of positive inte- 
gers {0,1,...}. 

2. For any £ > 1, a G S , the following decomposition holds 

I 



A.(0)^ + ^A.(rf(-a(z)|a) 



k=\ 



3. For any a G S, the following decomposition holds 



Ci(cj) = Mi 



^ oo 

A.(0)- + j;A.(A;)pf(-a(.)|a) 



fc=i 



(3.15) 



(3.16) 



Proof Item 1 follows directly from the definition and ()2.1 
To prove item 2, set 



cf(cT) = iM,A.(0) 



and observe that 



c^\a) = Y^[cf\a)-ct%)]+cf\a), 

k=l 

where the differences cf^^ {a) — cf^ (cr) are always positive. Since by definition 

r^\a)-ct'\<r)=MMk)I^\-<^{^)W), 
we obtain the decomposition stated in item 2. 

To prove item 3, observe that the uniform summability condition ()2.ip implies that 

lim cf^(o-) = Ci{cr) , 



(3.17) 
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\k] 

since — > for ^ oo. Thus, taking into account item 2, we obtain the desired 
decomposition. 



To present the coupling, it is convenient to rewrite the generator of the Glauber dynamics 
given in (I3.10p as 

^/(^) = EE^*(«l^)[/(^"")-/(^)]' (3.18) 

where 

= (^{j), for ah j / i, and cT*'''(i) = a, 

and where 

Ci{a) if a = —cr{i), 



c^{a\a) = \ ^_^^(^) ifa = a(.).' ^^'^^^ 



In a similar way, for any integer i > 1, we define the rates 

M.ELoA.(fc)-cP(^) ifa = cT(z 



and the generator 

^'"/(^) = E E^f'(«i^)[/(^^'") -/(^)]- (3-21) 

This amounts to include in the generators the rates of invisible jumps in which the spin i 
is updated with the same value it had before. Obviously the generator defined in (jS.lOp 
(and in (j3.12p ) defines the same stochastic dynamics as the generator given in (|3.18p (and 
in ()3.2ip . respectively). 

With this representation, we have the following corollary of Theorem [2j 



Corollary 1 



^/(^) = EEE^^*^*(^>f^(«i^)[/(^"")-/(^)]' 



and for any £ > 1, 



^"^/(^) = EEE^^*^^(^)?'l («i^)[/(^"")-/(^)]- 

ieZd a(^Ak=0 

From now on, we fix £ = L. The decomposition given in Corollary [T] suggests the follow- 
ing construction of the Glauber dynamics having Q and ^[^1 as infinitesimal generator. 
For each site i £ Z'^ consider a Poisson point process A^* having rate Mj. The Poisson 
processes corresponding to distinct sites are all independent. If, at time t, the Poisson 
clock associated to site i rings, we choose a range k with probability Aj(/c) independently 
of everything else. Then we update the value of the configuration at this site by choosing 
a symbol a with probability pf^\a\a{Bi{k))) depending only on the configurations inside 
the set Bi{k). It is clear that using this decomposition, we can construct both processes 
at and cjj'^' in a coupled way, starting from any initial configuration. 
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Actually we can do better than this. We can make a perfect simulation of the pair of 
measures and fi^^^ which arc the invariant probability measures of the processes having 

generators Q and Q^^^ respectively. 

Let us explain how we sample the configuration at a fixed site i G Z'^ under the measure /j,. 
The simulation procedure has two stages. In the first stage we determine the set of sites 
whose spins influence the spin at site i under equilibrium. We call this stage backward 
sketch procedure. It is done by climbing Tip from time back to the past a reverse time 
Poisson point process with rate Mj until the last time the Poisson clock rung. At that 
time, choose a range k with probability Xi{k). If A; = 0, we decide the value of the spin 
with probability ^ independently of anything else. If k is different from zero, we restart 
the above procedure from any of the sites j G Bi{k). The procedure stops whenever each 
site involved in this backward time evolution has chosen a range 0. 

When this occurs, we can start the second stage, in which we go back to the future 
assigning spins to all sites visited during the first stage. We call this procedure forward 
spin assignment procedure. This is done from the past to the future by using the update 
probabilities pf^ starting at the sites which ended the first procedure by choosing range 
0. For each one of these sites a spin is chosen by tossing a fair coin. The values obtained 
in this way enter successively in the choice of the values of the spins depending on a 
neighborhood of range greater or equal to 1. 

In the case of the measure /Lt'^l we use the same procedure, but only considering the choices 
of k which are smaller or equal to L. These procedures will be described formally in the 
next section. 

4 Perfect simulation of the measure /i and coupling 

4.1 Backward sketch procedure 

Before describing formally the two algorithms described at the end of the last section, let 
us define the stochastic process which is behind the backward sketch procedure. Our aim 
is to define, for every site i G Z*^, a process (Ci*^)s>o taking values in the set of finite 
subsets of Z'^, such that ci'-* is the set of sites at time —s whose spins affect the spin of 
site i at time t = 0. This is done as follows. 

For each i G Z'^, denote by . . . ^2 < T'i^ < < < < < . . . the occurrence times 

of the rate Mj Poisson point process iV* on the real line. The Poisson point processes 
associated to diff'crcnt sites arc independent. To each point Tj\ associate an independent 
mark iiT^ according to the probability distribution {Xi{k))k>o- As usual, we identify the 
Poisson point processes and the counting measures through the formula 

N'[s,t] = J2Ms<n<t}- 

The backward point process starting at time 0, associated to site i G Z'^ is defined as 

r^(^,o) ^ -Tl„+i, for any n>l. (4.22) 
We also define the associated marks 

= Ki^+, (4.23) 
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and for any k > 0, the backward A:-marked Poisson point process starting at time as 

N(^'%,u] = l^,<^(.o)<„^ V,.o)^,^. (4.24) 

ngZ 

To define the backward sketch process we need to introduce a family of transformations 
{7r*^*''^\ i G Z'^, /c > 0} on ^{11^), the set of finite subsets of Z'^, defined as fohows. For any 
unitary set {j} and A; > 1, 

For A; = 0, we define 

For any set finite set F dU^, we define similarly 

7r(^''=)(F) = U,-e^^(^''=)({j}). (4.27) 

The backward sketch process starting at site i at time will be denoted by (ci*^)s>o. The 
set Cs is the set of sites at time —s whose spins affect the spin of site i at time t = 0. 

The evolution of this process is defined through the following equation. Cq*'* = and 
/(CW) = /(C«) + r[fi^^'''\ci'l)) - fici'l)]N(^''^\du), (4.28) 

where / : T(T!^') ^ M is any bounded cylinder function. This family of equations charac- 
terizes completely the time evolution (ci*^)s>o. 

In a similar way, for the truncated process, we define its associated backward sketch process 
by 

/(Cfl'«) = /(Ci^l'«) + X; E /'[/(vr(^''=)(cfj'«)) -/(Cfj'«)]iv(^»(d^), (4.29) 

fc=0 

where we use the same Poisson point processes A^(J''=),0 <k<L, as in (j4.28p . 

The following proposition summarizes the properties of the family of processes defined 
above. 

Proposition 2 For any site i € Z'^, ci'^ and cl^^'^^\ s € M+, are Markov jump processes 
having as infinitesimal generator, 

^fiC) = Ej;M,A,(fc)[/(CUi?,(A:))-/(C)] + A,(0)[/(C\{i})-/(C)], (4.30) 

ieC k>l 

with initial condition at time t = 0, Cq^ = {i}, and 
L 

^"-'/(C) = 5]5]M,A.(A;)[/(CUi?i(A;))-/(C)] + Ai(0)[/(C7\{i})-/(C)], (4.31) 

i&C k=l 

with initial condition at time t = 0, Cq^''*'*'' = {i}, respectively. Here f : J^(Z'^) ^ R is 
any bounded cylindrical function. 
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r«oP = inf : = 0} 



Let 

- STOP 

~ (i) 

be the relaxation time. We introduce the sequence of successive jump times Tn ,n > 1, 
of processes N^^'^^ whose jumps occur in (|4.28p . Let T^^^ = T^^*'*^-* and define successively 
for n > 2 

r(^) = inf{t > : 3j G C% , 3k : N^^''^^ Qfj^l^tD = 1}. (4.32) 

Now we put 

C» = C«. (4.33) 

^ n 



Finally, notice that N^j^^Qp of Definition [41 is equal to 

7V«op = inf{n:C« = 0}. 
This is the number of steps of the backward sketch process. For the algorithm to be 

(i) 

successful, it is crucial to show that both relaxation time T^r^Qp and the number of steps 
^STOP finite. This is the content of next theorem. 

Proposition 3 Under the requirement 

sup^ \Bi{k)\Xi{k) < 1, (4.34) 

we have uniformly for any z G Z'^, 

PiT^TOP >t)< (4.35) 

and 

7 



E{N^lop) < 7,^ (4-36) 



where 

-f = l-snpJ2\Bi{k)\Xi{k). (4.37) 



The proof of Proposition [3] will be given at the end of Section [5] below. 
4.2 Algorithms 

We now go to the crucial point of precisely defining the maps F and G whose existence 
is claimed in Theorem [TJ These maps are implicitly defined by the backward sketch 
procedure and the forward spin assignment procedure. Let us give the algorithmic like 
description of these procedures. The following variables will be used. 

• is an auxiliary variables taking values in the set of non- negative integers {0, 1,2,...} 

• Ng}pQp is a counter taking values in the set of non-negative integers {0, 1,2,...} 

• Ng^Qp is a counter taking values in the set of non-negative integers {0, 1, 2, . . .} 
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• / is a variable taking values in 

• /I^^ is a variable taking values in Z*^ 

• if is a variable taking values in {0, 1, . . .} 

is a variable taking values in {0, 1, . . .} 

• S is an array of elements of Z*^ x {0, 1, . . .} 

is an array of elements of Z'^ x {0, 1, . . .} 

• C is a variable taking values in the set of finite subsets of Z*^ 

is a variable taking values in the set of finite subsets of 1^^ 

• is an auxiliary variable taking values in A 

• X is a function from Z'^ to A U {A}, where A is some extra symbol that does not 
belong to A 

is a function from Z"^ to A U {A}, where A is some extra symbol that does not 
belong to A 

We will present the backward sketch procedure only for G, the coupled perfect simulation 
algorithm. The backward sketch procedure for F can be immediately deduced ignoring 
the steps from 11. to 18. concerning the variables with superscript L. 

Algorithm 1 Backward sketch procedure 

1. Invut: {i}; Outvut: N^Iqp^ ^st'op' ^' ^^^^ 

2. N^O, N^op ^ 0. ^STOP ^ 0, i? ^ 0, ijl^l ^ 0, C ^ {i}, Ct^l ^ {i} 

3. WHILE C 7^ 

4. N^N+l 

5. Choose independent random variables ^jy' having exponential distribution of pa- 
rameter MjXj{k) for any j G C and for any A; > 

6. Choose {I,K) = aigmm{S^^'''\ j €C,k>0} 

7. IF K = 0,C ^C\{I} 

8. ELSE C ^ C U Bi{K) 

9. ENDIF 

10. B{N) ^ {I,K) 

11. IF 7^ 

12. ^ N 



12 



13. Choose = argmin{5S^'^\ j eC,0<k<L} 

14. IF = 0, Cl^l ^ \ {/W} 

15. ELSE Ct^l ^ Ct^l U (i^W) 

16. ENDIF 

17. bW(a^) ^ 

18. ENDIF 

19. ENDWHILE 

20. N^lop^N 

21. RETURN N^lop^ ^stop^ ^^^^ 

U B = SW, then we use the following Forward spin assignment procedure to sample 
X(i) = X^^\i). This is the algorithmic translation of the ideas presented in the last para- 
graph of Section 4. li B ^ B^^\ then we use the Forward spin assignment procedure twice 
independently, starting with input N^^Qp, B in order to sample X{i) and starting with 
input N^sToh -S'^' in order to sample 

Algorithm 2 Forward spin assignment procedure 

1. Input: N^sTOP' ^'^ Output: {X{i)} 

2. AT ^ iV« OP 

3. X{j) ^ A for all j G Z'^ 

4. WHILE AT > 1 

5. iI,K)^BiN). 

6. IF K = choose W randomly in A according to the probability distribution 

P{W = v)=pf\v) = ^ 

7. ELSE choose W randomly in A according to the probability distribution 



P{W = v) =pf\v\X) 



8. ENDIF 

9. X{I) ^ W 

10. N ^ N-l 

11. ENDWHILE 

12. RETURN {X{i)} 
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Remark 3 At a first look to steps 5 and 6 of Algorithm 1, the reader might think that the 
simulation of an infinite number of exponential variables is necessary in order to perform 
the algorithm. Actually, it is sufficient to simulate a finite number of finite valued random 
variables, see Knuth and Yao (1976) and also Section 10 of Galves et al. (2008). 



5 Proof of Theorem [T] 

This section is devoted to the proof of Theorem [TJ For the convenience of the reader, we 
start by recalhng the following theorem of Galves et al. (2008). 

Theorem 3 Under the requirement (I4.34p . 

sup V \Bi{k)\X^{k) < 1, 

the Gibbs measure fi is the unique invariant probability measure of the process {(Jt)t- The 
output X{i) obtained using successively Algorithms 1 and 2 given in Section\^ is a perfect 
sampling of the Gibbs measure fi. Finally, the estimate ^.35 ) holds : 



The same result holds true also for Xl^l obtained using successively Algorithms 1 and 2 
but now with restriction that only ranges k < L are considered. 

Remark 4 In Galves et al. (2008) the above result was proved for the minimal decom- 
position of the spin flip rates. The same proof works for the decomposition of the spin 
flip rates Cj(cj) considered in Section 4- Note that in contrary to the decomposition con- 
sidered in Galves et al. (2008), the decomposition in Section 4 is explicit in terms of the 
interaction. 

We first check that the assumption ()2.6p of Theorem [T] implies condition (j4.34p of Theorem 

El 

Lemma 1 Under the conditions of TheoremUl condition i4-34^ is satisfied. 
Proof For any A; > 0, write for short 

B:ieB,BglBi{k) 
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>fe-i 



For each i G Z'^ we have that 

oo oo 

Y,K{km{k)\<2dx,{i)+Y,mk) 

k=2 



k=l 



>k 



k=2 



-13 j: 



1) \Jb\ 



(5.38) 



B:i&B,B<ZBi{k),B<tBi{k-l) 



, fc=2 



Under the conditions of Theorem [H this last expression is finite. Therefore condition 
(I4.34P reads as 



B:ieB,BCBj^{l 



(oo 
fc=2 

Finally, let Pc be the solution of 

2dsup (e-^^»^' (l - e-^^^^-f'^^-^^'^W) 

i ^ ^ 

/ CO 

This concludes the proof. 



E i^b||<i. 

B:ieB,BcBi{k),B(^Bi{k-l) 



,k=2 



E I-^bI I = 1- (5.39) 

B:ieB,BcBi{k),B(lBdk-l) 



Proof of Theorem [H 

Proof of item 1. Theorem [3] combined with Lemma [T] proves that the output of the 
perfect sampling algorithm is indeed a sampling of /i. Next we prove (|2.7p . Define 

the cardinal of the set Cn^ after n steps of the algorithm. Let {Dn^)^^^ ^^^d be i.i.d. 
random variables, independent of the process, taking values in {—1,0, 1,2,...} such that 
for k > 1 

P{D^^ = \Bi{k)\-l) = K{k), (5.40) 

and for A; = 

P{D^^ = -1) = A,(0). (5.41) 

Note that under assumption (j2.6p . the condition (j4.34p holds (see Lemma 1). Hence, 
supjg^d i?(I?f ^) < —7 < 0. For any i G Z"^, let us call the index of the site whose 
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Poisson clock rings at the n—th jump of the process CW. Recah that /i'^ are conditionally 
independent given the sequence (Cn^)^ such that 



Put 



Note that by construction, 5n + n7 is a super-martingale. Then a very rough upper 
bound is 

L^i'^ < 1 + as long as n < V^tj^^p, (5.42) 
where V^^Qp is defined as 

F»OP = : sf = -1}. 

By construction 

^^STOP - ^STOP- 

Fix a truncation level N. Then by the stopping rule for super-martingales, we have that 

)+7i?(vi?opAiV)<0. 

But notice that 

E{S% ) = -1 • PiV^^p <N) + EiS^;f;V^%p > N). 

On {V^^op > N}, 5}^^ > 0, hence we have that E{S%^ ) > -P{V^top ^ We 
conclude that 



Now, letting ^ 00, we get 



and therefore 

milop) < ^- 

Proof of item 2. We have to show that the coupled perfect sampling algorithm G defined 
in Section 4 achieves the bound (j2.8p . For that sake we introduce the following stopping 
times. For any site j G Z"' and any t > 0, let 

Ti^^'*^ = inf{Ti^'0) >t}-t 
be the first jump of the Poisson point process (T^^'°^)„ shifted by time t. 



16 



Let us call the first time that a range of order k > L has been chosen. is given by 
t|*) = inf{i > : ^ ^ N^^'''\[t, t + t[^''\) > 1}. (5.43) 

Recall that in order to construct fi and we use the same Poisson point processes for 
all ranges k < L. Thus we have that 

P{X{i) + < P(T« < T«op)- (5.44) 

By Lemma [2] given below we conclude the proof of Theorem [1] • 



Lemma 2 



p(r« < r«op) < sup ( 



Proof Put 



a. 



{k) = Y^\{i). 



i<k 



Given the structure of the interaction, we have for k > 2, a-iik) = e ^^B:ieB.B(rBi{k)\-^B\ 
and thus 1 - ai{L) < sup^g^d (l - e"^^s:ies,fl!Zfl,(L) l-^sl^j _ CaU this last quantity 



Then we have 

p/rpii) ^ rp{i) \ 



6{L) = sup (l - e"^^s:>es,s!zs,(L) \Jb\\ _ 



<^e[p ({iV(^''^)(]r«i, r«]) = 1 for , G C« „ k>L}\ C« „ iviop > n - 1 



n>l 



M,(l-a,(L)) 



1, 



n>l 



< 



Mj6{L) 



n>l 



= 6iL)E{Nllop)- 

In the above calculus we used that given C^*l^, Tn^ defined in ()4.32p is a jump of N^^'^^ 
with probability 

MjXj{k) 

Summing over all possibilities k > L gives the term (Mj(l — aj(-Z>)))/(X^, M^). By 

(j4.36p we get the result. • 

Proof of Proposition [3] (j4.35p is the statement of Theorem [3l (|4.36p has been proved 
above in the proof of (|2.7p . • 
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6 Final discussion 



It is worth to note that Theorem [T] immediately provides an upper bound for the d— distance 
of fi and yul-^l. The J— distance is defined as fohows. 

Definition 5 Given two probability measures /i and u on S, a coupling between /x and v is 
a probability measure on S x S having as first and second marginals fi and v, respectively. 
The set of all couplings between /i and v will be denoted A4{fj,,i'). 

In the next definition, the elements of the product space S x S will be denoted by 
((ai(i),cT2(i)),i GZ'^). 

Definition 6 The distance d between two probability measures vi and V2 on S is defined 
as 

d{i^i,i^2) = inf < sup Q{ai{i) / fT2(i)) > . 

Q£M{yi,y2) [jg^d J 

This definition naturally extends Ornstein's d— distance to the space of non-homogeneous 
random fields. As a corollary of Theorem [T] we obtain 

Theorem 4 Under the assumptions of TheoremUl we have 
where C is the same constant as in (E 



Upper bounds on the d— distance can as well be obtained by using a contraction argument. 
This approach was introduced in Dobrushin (1968-1970) to prove uniqueness of the infinite 
volume Gibbs measure. Applied to our context, it implies the existence of a good coupling 
through compactness arguments. 

For a nice and self contained presentation of this method, we refer the reader to Pre- 
sutti (2009), see also Chapter 8 of Georgii. Applying Dobrushin's method, we obtain the 
following upper bound on Ornstein's d— distance. 

Theorem 5 Assume that 

/3sup ^ \Jb\ = r < 1. (6.46) 

Then there exist unique infinite volume Gibbs measures fi and fi^^^ and they satisfy 



^_ sup Yl (6.47) 

^^^'^ B:ieB,B<^B,(L) 



Proof The proof of Theorem can be achieved following the ideas in Dobrushin (1968- 
1970), exposed in Chapter 3 of Presutti (2009). • 
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Observe that condition (j6.46p is weaker than the one requested in Theorem[TJ In particular, 
this condition is too weak to guarantee that our perfect simulation procedure of /x and /i'^l 
stops after a finite number of steps. As a counterpart of the less restrictive assumption 
(j6.46p . the upper bound (j6.47p is not as good as the upper bound (|6.45p . 

However, from our point of view, the important fact is that the contraction method of 
Dobrushin is not constructive and does not provide an expHcit sampling procedure for both 
measures. To provide an explicit coupled perfect simulation procedure for the measures 
is precisely the goal of the present paper. 

Further features of high temperature region can be deduced by our constructive approach. 
We discuss as an example the decay of correlations between localized observables as the 
regions of localization are moved away from each other. 

This decay of correlations can be described in terms of an associated random walk that 
we introduce now. Suppose that is translational invariant. In this case, the sequence 
{Kik))k does not depend on i, and we shall write X{k) = Xi{k), for any k > 0. Let ^2, • • • 
be i.i.d. random variables with common distribution 

P(ei = -1)=A(0), P{^, = \Bo{k)\-l)=X{k), k>l. 

Put 93(A) = E{e^^^) and let 

p = sup{A > : ip{X) < 1}. 

Finally, consider the following left-continuous random walk Sn = ■Ci + • • • + Cn; for n > 0. 
Note that Sn is an upper bound on the total number of sites whose spins affect the 
spin at site i under equilibrium, after n steps of the algorithm. Under condition (j4.34p . 
-E(^i) = —7 < 0, hence Sn drifts to —00. This implies that M = sup„ Sn, the maximum 
of the random walk, is finite almost surely. 

The following theorem is an immediate consequence of our construction and of results 
of Korshunov (1997) concerning the maximum of a random walk. In order to state the 
theorem, we have to recall the notion of locally power functions. 

Definition 7 A function f : M4. — > IR+ is called locally power if for every fixed t > 0, 
f{x + t) ~ f{x) as X ^ 00, i.e. lima;^oo = 1- 

Theorem 6 Take two cylinder functions f and g having support Aj (s Z'^ and d 7/. 
Assume the distance (i(Aj, A^) = R. Under the assumption \4.34^ we have 

Hfg) - Kf)K9)\ < 2c/c,(|Ay| + 1A,|) P{M > R/2), (6.48) 

where 

c/ = supsup|/(a)-/(0| (6.49) 
and Cg is similarly defined. Concerning the decay of P{M > R/2), the following holds. 
1. Suppose that p = and that P{£,i + (,2 ^ n) 2P(,^i > n) as n ^ 00. Then as 

77, — > CX3, 

P{M > \Bo{n)\) ^ - 
7 

where 7 is given in (j4.37p . 



J2 X{k)\Bo{k)\ - \Bo{n)\ i X{k) 



_k=n+l 



\k>n+l 
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2. Suppose that p > 0, ip{p) < 1, that P{6.i + S,2 ^ n) ^ cP{S,i > n) as n ^ oo for some 
constant c and suppose that e^P{S,i > n) is locally power. Then as n ^ oo 

P(M>|5o(n)I)^f^— ( Y: ^(^)- 

k=n+l 

3. Suppose that p > 0, (p{p) = 1 and (p'{p) < +oo. Then 

P{M > \Bo{n)\) ~ ce-P^ 
as n —> oo, where c is some positive constant. 

Remark 5 Using the same type of upper bounds as in the proof of Lemma d note that 
the expression 

oo 

E ^(^)l^o(A:)| 

k=n+l 

that defines the decay of correlations can be upper bounded by 

k=n+l \B:OeB,BcBo{k),Bi^Boik-l) I 



which yields the tail of the converging series of condition 112. 6\) . 



Proof We start by proving (16.480 . It is sufficient to prove it for functions / and g having 
support = {i} and Ag = {j}, R = \\i — j\\ > 0. For any site k € Z*^, let 

n>0 



C« = IJ (6.50) 



be the total number of sites whose spins affect the spin of site k under equilibrium (recall 
(j4.33p ). Then it is evident that 

Hfg) - KfM9)\ < cfCg{\Af\ + |A,|)P(c« n + 0). (6.51) 

Then the same comparison argument that leads to (15.420 yields 

P(CW n C^-'') 7^ 0) < P(M + M' >R) < 2P{M > R/2). 
This concludes the proof of (|6.48p . 

Item 1 is a consequence of Theorem 1 of Korshunov (1997), see |10] . First notice that 
F{t) = P{(,i > t) is constant on intervals [|i?o(^)| ~ 1; |-So(n + 1)| — 1[. Then evidently 
F{\BQ{n) \ — 1) = 1 — a(n), where a„ = Aq + • • • + A^- Therefore, 

fOO 



/•oo / \ 

Fit)dt = [Y^(l-a{k))[\Bo{k + l)\-\Bo{k)\]] -{l-a{n)) 

= ( f; X{k)\Bo{k)\] - {1 - ainmBoin)\ + 1] 



\k=n+l / 

' oo \ / N 

Kk)\B,{k)\\-[\Bo{n)\ + l\[Ym 

\k=n+l J \k>n / 
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This integral F{t)dt is the quantity that determines the asymptotic behavior of 

P{M > \Bo{n)\) as it is shown in Theorem 1 of Korshunov (1997). Items 2 and 3 are a 
consequence of Theorem 2 and Theorem 3 (ii) of Korshunov (1997). • 
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